Efficient spin control in high-quality-factor planar micro-cavities 
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A semiconductor microcavity embedding donor impurities and excited by a laser field is modelled. 
By including general decay and dephasing processes, and in particular cavity photon leakage, detailed 
simulations show that control over the spin dynamics is significally enhanced in high-quality-factor 
cavities, in which case picosecond laser pulses may produce spin- flip with high-fidelity final states. 



I. INTRODUCTION 



Photons and excitons (X) can be made to strongly in- 
teract in high-quality cavities containing a semiconductor 
quantum well, leading to a repetitive coherent exchange 
of energy between the two particles J When the energy 
exchange ocurrs faster than the decay time of the individ- 
ual components, a combined state, the exciton-polariton, 
is said to have formed. 

Exciton-polaritons show a variety of features, that 
motivate studies in multiple directions. Current inter- 
est in exciton-polariton research in 2D micro-cavities 
mainly focus on its liquid state and non-equilibrium Bose- 
Einstein condensation (BEC) covering many aspects of 
the problem 2 . The interaction of polariton fields with 
impurities or defects has been studied for the case of a 
polariton fluid scattered by centers acting on the pho- 
tonic component of the field^, a polariton gas scattered 
by spin-independent disorder potential acting on the ex- 
citon degree of freedom^, and the scattering of polaritons 
from spinless impurities acting on the excitonic compo- 
nent of the fielc£. However, to the best of our knowledge, 
there are no reported studies on the interaction of a po- 
lariton field with a single spin degree of freedom. 

Here we study the dynamics, including relaxation pro- 
cesses, of a diluted exciton/photon field interacting with 
a single impurity of spin s = 1/2. The system is de- 
picted in Fig. [TJ It consists of a 2D photon cavity em- 
bedding a quantum well (QW), which contains few donor 
impurities^. The whole system is assumed to be at low 
temperature and excited by a laser from outside. We 
show that the quantum control of a single spin is more 
efficient for high-quality- factor cavities. Thus, a spin- 
flip in a high-fidelity final state could be produced with 
a single laser pulse of a few picoseconds. Since typi- 
cal decoherence times for impurity spins in semiconduc- 
tors are in the /jls time-scale, the system can act as a 
high speed quantum memory or qubit£— . We believe 
that the present proposal and that for the implementa- 
tion of two-qubit polariton-induced operation o 12 i 13 sug- 
gest that a complete quantum-computing scalable archi- 
tecture based on a solid-state system is possible using 
polaritons in 2D-microcavities. 
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FIG. 1. Pictorial representation of the system. Two Dis- 
tributed Bragg Mirror (DBR) structures, placed at the sides 
of a quantum well (QW), confine photons injected from out- 
side by a laser. The photons produce QW excitons that in- 
teract with the impurity spin localized at position R. 



II. COUPLING WITH THE ENVIRONMENT 

Real systems cannot be entirely isolated from their en- 
vironment. This is specially true for solid-state systems 
where several particles co-exist. When control is in mind, 
undesired interactions make the evolution unpredictable 
with the possibility of partial or total failure of the con- 
trol operation. In our case, excitons, photons, and spins 
suffer from the coupling with the environment. 

For highly pure samples, low temperatures and low ex- 
citon density, the relevant decoherence processes for ex- 
citons are those causing spin flip of electrons and holes, 
with conversion between bright and dark excito mj 18 i 22 . 
In general hole spins loss coherence faster than elec- 
tron spins; for instance, in CdTe the spin relaxation 
of electrons is 29ps^, while that for holes is < 7ps^i. 
In addition, the annihilation of excitons must be also 
considered, with associated lifetime of hundred ps in 
GalnNAs/GaAs^. 

While different processes, such as structural disorder- 
are responsible for the loss of photon population and co- 
herence, the main process is photon leakage off the cavity 
due to its finite Q-factor, which leads to a lifetime of the 
order of r = 15ps^. 
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At extremely low concentration of impurities with den- 
sities rii ~ 10 13 cm -3 , electrons bound to different 
donors are well localized and do not interact among 
them 27 . The interaction with the nuclei is dominant (due 
to the strong confinement of the localized state). At tem- 
peratures T < 10K the transverse relaxation time T2* 
is a few ns for the electron bound to a donor, and the 
spin relaxation time is of the order of /is for donors in 
GaAs 2 !. 



III. HAMILTONIAN 

In what follows, we work in the Heisenberg picture; 
thus, time-dependent operators shall be everywhere un- 
derstood. The free Hamiltonian reads 



k C x k C xk 



(i) 



where the first and second terms correspond to exci- 
tons (^k/^ak) and cavity photons (c %k /c x k). The QW 
quantum confinement splits the heavy- and light-hole 
electronic bands, forming excitons out of conduction- 
band electrons with total angular momentum j z = 1/2 
and valence-band heavy holes with j z = 3/2. Bright 
(j z = 1) and dark (j z = 2) excitons are included: a = 
{1,2,3,4} = {tl^lt^4} 5 where the sin S le (double) 
arrow identifies an electron (hole) angular momentum. 
The respective dispersion relations for excitons and pho- 
tons are = £o + (hk) 2 /2m* and = c/n (k 2 ^k 2 ) 1 / 2 , 
where k is the in-plane momentum; the momentum k 2 
in the growth direction is determined by parameters of 
the cavity, n is the index of refraction, and c the speed of 
light. The polarization of the photon is x — {1, 2}. The 
ground state energy of the donor is set to zero. 

The system is excited by a classical laser field pro- 
ducing photons that propagate inside the cavity. Using 
the quasi-mode approximation (useful for high Q-factor 
cavities)^, the cavity-laser interaction reads 

H LC = hVI^ V xk (t) e^*-^ c xk + H.c. , (2) 

Xk 

where V xq (£) is the coupling constant, A = L 2 is the area 
of the system, is the laser frequency and Cl a constant 
adequately chosen to ease the numerical solution, see be- 
low. 

Cavity photons interact with excitons according to 

H L = ^2 ^xk(^k)c xk ^ k + i^.c. , (3) 

X«k 

where ^ Q!X k(^k) = for a = 1,4. 

The QW contains donor impuritites. The assumption 
is made that, at low temperature, each impurity has an 
electron bound to it that contributes a spin s = 1/2; in 
addition, the concentration of donors is low enough to 
ensure that excitons/polaritons will interact only with 



one selected impurity, located at position R, when the 
laser spot is small enough* 2 ^ Via Coulomb exchange, the 
electrons belonging to the exciton and the donor impurity 
interact through Hxs = Hx~s + H- c - 
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where the vector spin operator s = (s x ,s y ,s z ), jkk' = 
jo[l + aj 2 (k — k 7 ) 2 ] -1 / 2 , and ai is a measure of the im- 
purity electron localization^. Here we adopted a more 
detailed notation for the spin a of the exciton: the elec- 
tron (hole) spin has index x (v)i an d cr is the vector of 
Pauli matrices. 9 X-X interaction is disregarded, because 
we will study the case of low exciton concentration, where 
nxa*B 1 / ^ < 1, with the exciton Bohr radius. 



IV. METHOD 

Different theoretical tools are employed to solve prob- 
lems in exciton-polariton research. Heisenberg equa- 
tions of motion describe the dynamics of mean val- 
ues of either exciton/photon operators or polariton 
operators 3 - ! 30 ! 31 . It is also common the use of the Gross- 
Pitaevskii equatio n 25 ! 26 . Other methods have also been 
used, such as the Hartree-Fock-Popov 3 - 2 -. 

We make use of the Heisenberg equations of motion 



(HEM) hd(0)/dt = i( 



H,0 



) for mean values ((...)) of 



operators describing separately excitons, photons, and 
the impurity spin. This allows us to treat the cases of 
weak — where no polaritons exist — and strong coupling, 
as well as to include easily spin-flip processes that cause a 
polariton to dissociate into a dark exciton and a photon. 

In general, the HEM comprise a set of infinitely 
coupled equations, that can be ordered in a heirar- 
chy, much as the Bogoliubov-Born- Green-Kir kwood- 
Yvon (BBGKY) hierarchy of classical statistical mechan- 
ics. One must then set equations for products or cor- 
relation of increasing number of operators. In order 
to close the system of equations, a truncation of the 
hierarchy is necessary. We use the truncation scheme 
(Oi0 2 ) = (d 1 )(d 2 }. It is important to note that the 
photon-exciton coupling is not affected by the truncation 
scheme; therefore, the formation of polaritons (strong- 
coupling regime) is accurately described. 

The system-bath coupling is properly introduced in the 
HEM by the formalism of Quantum Heisenberg-Langevin 
equations, which leads to additional terms in the equa- 
tions: damping, Lamb shift, and stochastic force T . A 
simpler way to deal with the environment is by introduc- 
ing constants, taken from experiments or other theoreti- 
cal works, directly in the HEM, taking into account the 
results of the detailed microscopic derivation. Here we 
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follow the phenomenological procedure, by adding con- 
stants directly in the HEM^ 3 -. Photons are coupled to 
the radiation field at zero temperature outside of the cav- 
ity, resulting in the addition of a term — £ q c xq ; the damp- 
ing £ q becomes very large when q is such that the normal 
component of the field exceeds the critical angle sepa- 
rating low and high DBR reflectivity. For excitons we 
introduced a term —/3 a b a ^ with a spin-dependent con- 
stant f3 a , accounting for radiative recombination and spin 
flip (no scattering is considered). For the impurity spin, 
a general constant 7 is used in all component, since no 
extenal magnetic field exists to distinguish among them. 
For long times, the spin relaxes, but does not vanish; 
thus, an equilibrium state is defined. 

We consider a circularly-polarized laser field that ex- 
cites 620 5 an d an impurity located at R = 0. To simplify 
the calculations, we eliminate fast oscillations by moving 
to a rotating reference frame, with frequency O, setting 
(6) = e~ int (&}, for O = c,b. For the sake of simplic- 
ity, we hereafter denote the rotating frame version (0'), 
simply as O. The equations of motion read 

ds z 
~dt 

ds x 
~dt 

dSy 

~dt 

where s z = 
n x 

rnk,mk' 
^(Pnkjmk' 
<?kk' = ^kk' 



kk' 

~A ^kk' ( S 2/Plk, 



5 ^lk,2k 
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- Pmk,nk')/ 2 r^ with p nk ,mk' = ^* k &mk' and 
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where s± = s x ± is y . Similar equations hold between 
&2q fe q and 6i q 6 4q . 
dr 
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V. RESULTS 

Numerical solution of the HEM is obtained using a 4th- 
order Runge-Kutta method in a 2D grid of TV x TV modes 



in momentum space. Basic units are {meV, ps, nm}, and 
we use data compatible with GaAs 35 . The values of the 
different parameters are taken, in most cases, directly 
from experimental or theoretical work, only Vo and jo 
are adjusted using our calculations. When g a ' a<l = 0, 
^2q(0) 7^ and s z (0) = ft, the system of equations be- 
comes linear, and can be solved exactly, jo is then ad- 
justed to yield a negative eigenvalue that matches the 
reported binding energy of excitons to donors (about 
lmeV). The value so obtained for h 2 jo/A ~ 10 -5 meV is 
in agreement with previous reports^^. We fix the value 
of the coupling Vo by demanding that the total exciton 
density n x = ^ q &J q &*q be low, i. e.: r = n x a*£/A < 1, 
so that the X-X interaction can be neglected. 

We studied the evolution of spin components, exci- 
ton and photon populations when the system, repre- 
sented by TV = 50 modes (larger TVs do not change 
the result significantly), is excited by a circularly polar- 
ized normal-incidence monochromatic laser-pulse Vcrk = 
Vo exp{ — (t — t p ) 2 /w 2 } . Cases with and without decoher- 
ence are considered. 

It is instructive to analyze first the (idealized) 
decoherence-free case — no plot presented. We find nei- 
ther &3 q 63 q nor &4 q &4q populations, while s z and fri q fri q 
change little from their initial values. The small change 
in s z can be understood as follows: according to Eqs. 
(O ds z /dt oc (ftjo/A)6^ q 6 2q and to Eqs. flS) dbi^/dt oc 
(ftj /i4)&2q, that roughly yields ds z /dt oc (ftjoM) 2 &2 q &2q- 
This, compared to ds x /dt oc (7yo/^)&2q&2q> is a very 
small quantity given the choosen value of hjo/A ~ 10 -5 
ps -1 . On the contrary, the spin projection in the xy- 
plane can rotate several cycles depending on the tempo- 
ral width and intensity of the pulse. As it is well known 
from quantum optics, once the laser is turned off, there is 
a remaining oscillating population of excitons and pho- 
tons. For certain values of the pulse parameters, these 
populations are so small (r —> 0) that cannot produce 
important changes in the spin. 

When decoherence is included, there is conversion to 
dark states &4 q , due to hole spin-flip, and to a lesser 
extent due to electron spin-flip. Because of the long 
life-time of these dark states, the fraction r remains fi- 
nite (though very small compared to its peak value). 
Fig. [2] presents the results for a simulation with pa- 
rameters {^0 = 2270.ps- 1 , Q = 2301.2ps" 1 ,e /ft = 
2301.2ps- 1 ,£ = 6.610- 2 ps-\/?i = /3 4 = 0, fa = £3 = 
lO^ps- 1 ,/^ = f3 34 = 3W- 2 ptT 1 ,p 13 = faA = lps- 1 }. 
We find that if Vo < 32ps _1 nm _1 then r < 1 and the ne- 
glect of the X-X interaction is justified. Under this con- 
dition, we see that a single inversion s x — >> — s x can be re- 
alized in few picoseconds. Faster spin motion is observed 
when the laser intensity (and so the photon/exciton pop- 
ulations) increases. As predicted in the previous para- 
graph, during the whole evolution, the change in s z is 
very small, as seen in the lower panel of Fig. [2j In ad- 
dition the population of dark excitons is also very small 
compared to that of bright excitons: with the definition 
r i = b*^iqCL*B/A, we obtain at t = lOps {r\ ~ 0.3, r2 — 
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1.5 x 10- 6 } and at t = 20ps {n ~ 6 x lO" 6 , r 2 ~ lO" 6 }. 
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FIG. 2. (color online) Evolution under the excitation by a 
laser pulse of width w ~ 4.5ps and Vo = 25ps -1 nm -1 ; the 
initial state is a spin having mean values {s x = y/3h/A, s y = 
0, s z — —ft/4}. Upper Panel: Spin components s x (solid blue) 
and s y (red). Lower Panel: Spin component s z . Inset: fraction 



A. Spin rotation for strong and weak coupling 

The addition of decoherence allows us to address the 
regimes of strong and weak coupling, and in particular to 
study the effect that cavity losses have in the spin con- 
trol. Weak coupling is characterized by |£o — fii \ > 2g/h 
(in our case 2g/h ~ 2.2 ps _1 ), and this regime can be 
simulated by increasing the photon losses of all modes 
(increasing £o)> which amounts to considering different 
cavities with varying quality factor Q. Two notes of cau- 
tion: First, we have treated the laser-photon coupling 
in the quasimode approximation, valid for high-quality- 
factor cavities. Therefore, we will refrain from studying 
cases with large values of £o- Second, the laser-photon 
coupling Vo is, in general, affected by changes in the pho- 
ton losses £o; however, we can envisage situations where 
one can increase £o without affecting Vo, for example - 
but not exclusively- by reducing only the reflectivity of 
the left DBR in Fig. [Q 

Fig [3] shows the effect that the increase in photon leak- 
age, at fixed laser field intensity, has on the rotation of 
the impurity spin. For simplicity other sources of deco- 
herence are disregarded. For all simulations, we used one 



set of laser parameters {w ~ 4.5ps, Vo = 15ps -1 nm -1 } 
for a gaussian pulse which produces, without photon loss, 
a rotation from the initial state s x = \/3ft/4 to the final 
state s x = — \/3ft/4 at t = 15ps, i. e. a change in the 
angle AO = tt. Next we simulated situations of increas- 
ing £o and plotted A0(£o)- In addition, we plotted the 
maximum photon population acchieved during the pulse. 
We observe that for £o < 2 ps -1 (Q > 1500) there is al- 
most full rotation of s x , and that for lower quality- factor 
cavities (high £o) the spin changes little. The population 
of cavity photons and excitons follow this tendency. 
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FIG. 3. (color online) Degree of spin in-plane rotation AO 
(dashed blue) as a function of the cavity photon loss, and 
maximum cavity-photon population (solid red). Inset: Zoom- 
in of rotation angle AO for low £o- 

The effect of decoherence can be characterized with the 
fidelity F. If the final state we wish to obtain is the pure 
spin state -1/2 |t) + y/3/2 ||), having mean values {s x — 
— V / 3^/4,s 2y = 0,s z = ft/4}, the formula for the fidelty 
reduces to F = (— l/2h)(V3s x — s z — ft), see Jozsa 36 . For 
the cases £o = 3.5ps _1 and £o = 20.5ps _1 the resulting 
fidelity is F = 0.9965 and F = 0.607, respectively. 

We interpret the enhanced rotation in high-Q cavities 
in the following way. For high Q, as seen in Fig. [3j the 
photon density is larger, and a repetitive and longer in- 
teraction with excitons is possible. This leads concomi- 
tantly to the formation of polaritons, with the excitonic 
component causing impurity spin rotations. In contrast, 
for lower values of Q, photons tend to leave the cavity 
faster, and there is small conversion to excitons. As was 
mentioned before, it is perhaps easier to envisage a cav- 
ity, whose Q factor is lowered by degrading the left DBR 
in Fig. [TJ Then, a naive picture tells us that the laser 
field produces photons inside the cavity at the same rate 
in the high and moderate Q cases. In the latter, photons 
are more prompt to leak out and produce less excitons. 

In addition, we can ask what the pulse width w should 
be to ensure a full rotation of s x , for different values of 
cavity loss (see Fig. 3]). As expected from the previous 
analysis, we see that one requires longer pulses to pro- 
duce the rotation, but in contrast to what happened be- 



5 



fore, the fidelity is almost unchanged. We attribute the 
behavior of F to the fact that the only source of decoher- 
ence is photon loss in these simulations and that the final 
state is forced (by changing w) to be the closest possible 
to the ideal state. 
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FIG. 4. (color online) Pulse width w (solid blue) required to 
produce full rotation of the s x spin component, and corre- 
sponding fidelity F (dashed red) as a function of cavity loss 

ft. 

Finally, it is worth mentioning that we have used a con- 
servative value for jo. For example, Puri et at^ reports 
a much higher value of jo for QDs replacing impurities. 
This would lead to even faster spin control, together with 
more efficient control on s z . However, fs laser pulses have 
a broad frequency spectrum, and may excite several po- 
lariton modes. This may lead to destructive interference 
effects, which may reduce the effectiveness of a large jo. 



VI. CONCLUSIONS 



We studied the optical control of single spins in micro- 
cavities accounting for all sources of decoherence. When 
the system is in the strong-coupling regime, the spin ma- 
nipulation is most efficient and can be done by a few- 
picoseconds laser pulse. This suggests that single spins 
embedded in high Q-factor planar cavities can act as 
quantum memories and as qubits, with the optical ex- 
citation being the mechanism to control the state of the 
memory or to perform one-qubit operations for quantum 
computing. This optical control produces high-fidelity fi- 
nal states in very short time: a single operation can be 
performed 10 6 faster than the typical decoherence time of 
the impurity spin qubit (compared to other proposals us- 
ing for example ion traps^I) with a fidelity of F > 99.8%. 
We believe that the present proposal for one-qubit opera- 
tions, together with a previous one for implementing two- 
qubit operations in the same syste m 12 ! 13 show that po- 
laritons in 2D-microcavities is a promising system for the 
implementation of solid-state quantum-computing scal- 
able architectures. 
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